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Abstract 



This paper is our third step towards developing a theory of testing 
^ ' monomials in multivariate polynomials and concentrates on two problems: 

00 ' (1) How to compute the coefficients of multilinear monomials; and (2) how 

t~^ , to find a maximum multilinear monomial when the input is a IIEn poly- 

^^D ' nomial. We first prove that the first problem is #P-hard and then devise a 

Cn I 0*(3"s(n)) upper bound for this problem for any polynomial represented 

r~^ . by an arithmetic circuit of size s{n). Later, this upper bound is improved 

f^ ' to 0*(2") for flEn polynomials. We then design fully polynomial-time 

("^ , randomized approximation schemes for this problem for IlE polynomials. 

On the negative side, we prove that, even for flEII polynomials with terms 
of degree < 2, the first problem cannot be approximated at all for any 
approximation factor > 1, nor "weakly approximated" in a much relaxed 
J^ ' setting, unless P=NP. For the second problem, we first give a polynomial 

?-i I time A-approximation algorithm for HEn polynomials with terms of de- 

grees no more a constant A > 2. On the inapproximability side, we give 
a n'^~^' " lower bound, for any e > 0, on the approximation factor for 
HEn polynomials. When terms in these polynomials are constrained to 
degrees < 2, we prove a 1.0476 lower bound, assuming P ^ NP; and a 
higher 1.0604 lower bound, assuming the Unique Games Conjecture. 

1 Introduction 

1.1 Background 

We begin with two examples to exhibit the motivation and necessity of the 
study about the monomial testing problem for multivariate polynomials. The 




first is about testing a fc-path in any given undirected graph G = {V, E) with 
\V\ = n, and the second is about the satisfiability problem. Throughout this 
paper, polynomials refer to those with multiple variables. 

For any fixed integer c > 1, for each vertex Vi £ V , define a polynomial pk.i 
as follows: 

Pi,i = 

Pk+l,i = Xl\ > ■ Pk,j I , fc > 1. 

We define a polynomial for G as 

n 
i=l 

Obviously, p{G^ k) can be represented by an arithmetic circuit. It is easy to see 
that the graph G has a fc-path Vi^- ■ ■ Vi^ iS p{G, k) has a monomial x'l ■ ■ ■ x1 of 
degree ck in its sum-product expansion. G has a Hamiltonian path iff p{G, n) 
has the monomial x^ ■ ■ ■ x'^ oi degree en in its sum-product expansion. One can 
also see that a path with some loop can be characterized by a monomial as well. 
Those observations show that testing monomials in polynomials is closely related 
to solving fc-path, Hamiltonian path and other problems about graphs. When 
c = 1, .Tij • • • ccjj. is multilinear. The problem of testing multilinear monomials 
has recently been exploited by Koutis [H] and Williams [3D] to design innovative 
randomized parameterized algorithms for the fc-path problem. 

Now, consider any CNF formula / ~ /i A ■ • • A /,„, a conjunction of m 
clauses with each clause fi being a disjunction of some variables or negated 
ones. We may view conjunction as multiplication and disjunction as addition, 
so / looks like a "polynomial", denoted by p{f). p{f) has a much simpler HE 
representation, as will be defined in the next section, than general arithmetic 
circuits. Each "monomial" tt = tti . . . tt™ in the sum-product expansion of p{f) 
has a literal tt^ from the clause fi. Notice that a boolean variable x e Z2 has 
two properties of a;^ = x and xx = 0. If we could realize these properties 
for p{f) without unfolding it into its sum-product, then p{f) would be a "real 
polynomial" with two characteristics: (1) If / is satisfiable then p{f) has a 
multilinear monomial, and (2) if / is not satisfiable then p{f) is identical to 
zero. These would give us two approaches towards testing the satisfiability of /. 
The first is to test multilinear monomials in p{f), while the second is to test the 
zero identity oi p{f). However, the task of realizing these two properties with 
some algebra to help transform / into a needed polynomial p{f) seems, if not 
impossible, not easy. Techniques like arithmetization in Shamir |28j may not be 
suitable in this situation. In many cases, we would like to move from Z2 to some 
larger algebra so that we can enjoy more freedom to use techniques that may not 
be available when the domain is too constrained. The algebraic approach within 
.^2(^1] in Koutis [21] and Williams [SOj is one example along the above line. It 



was proved in Bshouty et al. [B] that extensions of DNF formulas over ZJ to 
Zjv-DNF formulas over the ring Z^ are learnable by a randomized algorithm 
with equivalence queries, when A^ is large enough. This is possible because a 
larger domain may allow more room to utilize randomization. 

There has been a long history in theoretical computer science with heavy in- 
volvement of studies and applications of polynomials. Most notably, low degree 
polynomial testing/representing and polynomial identity testing have played in- 
valuable roles in many major breakthroughs in complexity theory. For example, 
low degree polynomial testing is involved in the proof of the PCP Theorem, 
the cornerstone of the theory of computational hardness of approximation and 
the culmination of a long line of research on IP and PCP (see, Arora at el. [3] 
and Feige et al. [HI). Polynomial identity testing has been extensively stud- 
ied due to its role in various aspects of theoretical computer science (see, for 
examples, Chen and Kao [H], Kabanets and Impagliazzo [T^) and its applica- 
tions in various fundamental results such as Shamir's IP=PSPACE [55] and the 
AKS Primality Testing [5]. Low degree polynomial representing [55] has been 
sought for so as to prove important results in circuit complexity, complexity 
class separation and subexponential time learning of boolean functions (see, for 
examples, Beigel |S|, Fu[TS], and Klivans and Servedio HO])- These are just a 
few examples. A survey of the related literature is certainly beyond the scope 
of this paper. 

1.2 The First Two Steps 

The above two examples of the /c-path testing and satisfiability problems, the 
rich literature about polynomial testing and many other observations have mo- 
tivated us to develop a new theory of testing monomials in polynomials repre- 
sented by arithmetic circuits or even simpler structures. The monomial testing 
problem is related to, and somehow complements with, the low degree testing 
and the identity testing of polynomials. We want to investigate various com- 
plexity aspects of the monomial testing problem and its variants with two folds 
of objectives. One is to understand how this problem relates to critical prob- 
lems in complexity, and if so to what extent. The other is to exploit possibilities 
of applying algebraic properties of polynomials to the study of those critical 
problems. 

As a first step towards testing monomials, Chen and Fu |5] have proved 
a series of results: The multilinear monomial testing problem for IIEn poly- 
nomials is NP-hard, even when each clause has at most three terms and each 
term has a degree at most 2. The testing problem for HI] polynomials is in 
P, and so is the testing for two-term nSII polynomials. However, the testing 
for a product of one two-term nSII polynomial and another IIS polynomial is 
NP-hard. This type of polynomial products is, more or less, related to the poly- 
nomial factorization problem. We have also proved that testing c-monomials 
for two-term IISII polynomials is NP-hard for any c > 2, but the same testing 
is in P for nS polynomials. Finally, two parameterized algorithms have been 
devised for three-term IIEII polynomials and products of two-term IIEII and 



nS polynomials. These results have laid a basis for further study about testing 
monomials. 

In our subsequent paper, Chen at al. [S] present two pairs of algorithms. 
First, we prove that there is a randomized 0*{p^) time algorithm for testing 
p-monomials in an n-variate polynomial of degree k represented by an arith- 
metic circuit, while a deterministic 0*{&A^ + P^) time algorithm is devised 
when the circuit is a formula, here p is a given prime number. Second, we 
present a deterministic 0*{2^) time algorithm for testing multilinear monomi- 
als in 11^22114 X Ilfclls polynomials, while a randomized 0*{l.b^) algorithm is 
given for these polynomials. The first algorithm extends the recent work by 
Koutis [H] and Williams [30] on testing multilinear monomials. Group alge- 
bra is exploited in the algorithm designs, in corporation with the randomized 
polynomial identity testing over a finite field by Agrawal and Biswas [1] , the de- 
terministic noncommunicative polynomial identity testing by Raz and Shpilka 
pS] and the perfect hashing functions by Chen at el. [TT]. Finally, we prove 
that testing some special types of multilinear monomial is W[l]-hard, giving 
evidence that testing for specific monomials is not fixed-parameter tractable. 

1.3 Contributions 

Naturally, testing for the existence of any given monomial in a polynomial can be 
carried out by computing the coefficient of that monomial in the sum-product 
expansion of the polynomial. A zero coefficient means that the monomial is 
not in the polynomial, while a nonzero coefficient implies that it is. Moreover, 
coefficients of monomials in a polynomial have their own implications and are 
closely related to central problems in complexity. As we shall exhibit later, the 
coefficients of multilinear monomials correspond to counting perfect matchings 
in a bipartite graph and to computing the permanent of a matrix. 

Consider a nSII polynomial F . F may not have a multilinear monomial in its 
sum-product expansion. However, one can always find a multilinear monomial 
via selecting terms from some clauses of F, unless all the terms in each clause 
of F are not multilinear or F is simply empty. Here, the real challenging is how 
to find a longest multilinear from the prod of a subset of clauses in F. This 
problem is closely related to the maximum independent set, MAX-k-2SAT and 
other important optimization problems in complexity. 

Because of the above characteristics of monomial coefficients, we concentrate 
on two problems in this paper: 

1. How to compute the coefficients of multilinear monomials in the sum- 
product expansion of a polynomial? 

2. How to find/ approximate a maximum multilinear monomial when the in- 
put is a HEH polynomial? 

For the first problem, we first prove that it is #P-hard and then devise 
a 0*(3"s(n)) time algorithm for this problem for any polynomial represented 
by an arithmetic circuit of size s{n). Later, this 0*(3"s(n)) upper bound is 



improved to 0*(2") for ZISII polynomials. Two easy corollaries are derived 
directly from this 0*(2") upper bound. One gives an upper bound that matches 
the best known 0*(2"') deterministic time upper bound, that was due to Ryser 
[26] in early 1963, for computing the permanent of an n x n matrix. The other 
gives an upper bound that matches the best known O* (1.415") deterministic 
time upper bound, that was also due to Ryser [5S], for counting the number of 
perfect matchings in the a bipartite graph 

We then design three fully polynomial-time randomized approximation schemes. 
The first approximates the coefficient of any given multilinear monomial in a HT, 
polynomial. The second approximates the sum of coefficients of all the multilin- 
ear monomials in a HE polynomial. The third finds an e-approximation to the 
coefficient of any given multilinear monomial in a IlfcEallt x Hm^s polynomial 
with a being a constant > 2. 

On the negative side, wc prove that, even for IISII polynomials with terms 
of degree < 2, the first problem cannot be approximated at all regardless of the 
approximation factor > 1. We then consider "weak approximation" in a much re- 
laxed setting, following our previous work on inapproximability about exemplar 
breakpoint distance and exemplar conserved interval distance of two genomes 
PUIIT]. We prove that, assuming P ^ NP, the first problem cannot be approxi- 
mated in polynomial time within any approximation factor a{n) > 1 along with 
any additive adjustment j3{n) > 0, where a{n) and /3(n) are polynomial time 
computable. 

For the second problem, we first present a polynomial time A-approximation 
algorithm for IIEn polynomials with terms of degrees no more a constant A > 
2. On the inapproximability side, we give a n^^~'^^" lower bound, for any 
e > 0, on the approximation factor for IIEII polynomials. When terms in these 
polynomials are constrained to degrees < 2, we prove a 1.0476 lower bound, 
assuming P ^ NP. We also prove a higher 1.0604 lower bound, assuming the 
Unique Games Conjecture. 

1.4 Organization 

The rest of the paper is organized as follows. In Section 2, we introduce the nec- 
essary notations and definitions. In Section 3, coefficients of multilinear mono- 
mials in polynomials arc shown to be related to perfect matchings in bipartite 
graphs and to the permancnts of matrices. Two parameterized algorithms are 
devised for computing the coefficient of a multilinear monomial with applica- 
tions to counting perfect matchings and computing the permanent of a matrix. 
In Section 4, we design three fully polynomial-time randomized approximation 
algorithms. Sections 5 and 6 are devoted to inapproximability and weak inap- 
proximability for computing multilinear monomial coefficients. Section 7 focuses 
on the problem of finding a maximum multilinear monomial in a polynomial. 
One approximation algorithm and three lower bounds on approximation factors 
are included. 



2 Notations and Definitions 

For variables xi, . . . ,Xn, let 'P[xi, • ■ • , x„] denote the communicative ring of all 
the n-variate polynomials with coefficients from a finite field V. For 1 < ii < 
■ ■ ■ < ik < n, TT = x^i^ ■ ■ ■ x]'' is called a monomial. The degree of tt, denoted 

by deg(7r), is X)s=iis- ^ i^ multilinear, if ji = • ■ • = j/t = 1, i.e., tt is linear 
in all its variables Xi^,. . . ,Xj^. For any given integer r > 1, tt is called a r- 
monomial, if 1 < ji, ■ ■ ■ ,jk < t. In the setting of the MAX-Multilinear Problem 
in Section [71 we need to consider the length of the a monomial tt ~ x-'^ ■ ■ ■ x-'-'' as 

\tt\ = X^fci log(l + jf). (Strictly speaking, |7r| should be J2e=i log(l + jf) logn. 
But, the common logn factor can be dropped for ease of analysis.) When tt is 
multilinear, |7r| = fc, i.e., the number of variables in it. 

For any polynomial F{xi, . . . ,Xn) and any monomial tt, we let c{F, tt) denote 
the coefficient of n in the sum-product of F, or in F for short. If tt is indeed in 
F, then c{tt) > 0. If not, then c{F,Tr) = 0. We also let S{F) denote the sum of 
the coefficients of all the multilinear monomials in F. When it is clear from the 
context, we use c(7r) to stand for c{F,tt). 

An arithmetic circuit, or circuit for short, is a direct acyclic graph with + 
gates of unbounded fan-in, x gates of fan-in two, and all terminals corresponding 
to variables. The size, denoted by s{n), of a circuit with n variables is the 
number of gates in it. A circuit is called a formula, if the fan-out of every gate 
is at most one, i.e., its underlying direct acyclic graph is a tree. 

By definition, any polynomial F{xi, . . . , Xn) can be expressed as a sum of a 
list of monomials, called the sum-product expansion. The degree of the polyno- 
mial is the largest degree of its monomials in the expansion. With this expres- 
sion, it is trivial to see whether F{xi, . . . ,a;„) has a multilinear monomial (or 
a monomial with any given pattern) along with its coefficient. Unfortunately, 
this expression is essentially problematic and infeasible to realize, because a 
polynomial may often have exponentially many monomials in its expansion. 

In general, a polynomial F{xi, . . . ,Xn) can be represented by a circuit or 
some even simpler structure as defined in the following. This type of repre- 
sentation is simple and compact and inay have a substantially smaller size, 
say, polynomially in n, in comparison with the number of all monomials in the 
sum-product expansion. The challenge is how to test whether F has a multilin- 
ear monomial, or some other needed monomial, efficiently without unfolding it 
into its sum-product expansion? The challenge applies to finding coefficients of 
monomials in F. 

Throughout this paper, the 0*{-) notation is used to suppress poly(ri, fc) 
factors in time complexity bounds. 

Definition 1 Let F{xi, . . . , Xn) S 'P[xi-, . . . , Xn] be any given polynomial. Let 
m,s,t > 1 be integers. 

• -F(a;i, . . . , a;„) is said to be a Tlm^s^t polynomial, if F{xi, . . . ,Xn) = 
Yii^iPij Fi — X^iLi^ij '^"'^ 1 < ''i < s, and Xij is a product of vari- 
ables with deg{Xij) < t. We call each Fi a clause. Note that Xij is not 



a monomial in the sum-product expansion of F{xi, . . . ,a;„) unless m = 1. 
To differentiate this subtlety, we call Xij a term. 

• In particular, we say F{xi, . . . ,x„) = Y[i=i ^i ^^ '^ ^m^s polynomial, if 
it is a HnjEsIIi polynomial. Here, each clause in F^ is a linear addition 
of single variables. In other word, each term in Fi has degree 1. 

• F[xi, . . . ,Xn) is called a nmSjIIt x IlfcE^ polynomial, if F{xi, . . . ,x„) = 
Fi ■ F2 such that Fi is a Hm^s^t polynomial and F2 is a Hk'^e polynomial. 

When no confusion arises from the context, we use IIEn and nS to stand 
for XlmSsIIt and IlmSs. respectively. 

Similarly, we use XIEsII and XlSg to stand for Hm'Ss^t and IlmSs respec- 
tively, emphasizing that every clause in a polynomial has at most s terms or is 
a linear addition of at most s single variables. 

3 Multilinear Monomial Coefficients, Perfect Match- 
ings and Permanent s 

In this section, we show that the problem of computing the coefficients of mul- 
tilinear monomials in a IISII polynomial is closely related to the problem of 
counting the number of perfect matchings in a bipartite graph and to the per- 
manent of a matrix with nonnegative entries. We first shall prove that com- 
puting the coefficient of any given multilinear monomial in a IISII polynomial 
is :/^P-hard. We then devise a 0*(3" s{n)) time fixed parameter algorithm for 
computing coefficients for multilinear monomials in a polynomial represented 
by an arithmetic circuit of size s{n). This upper bound is further improved to 
0*(2") for nEn polynomials. As two simply corollaries of this latter upper 
bound, we have an 0*(1.45") to find the number of perfect matchings in any 
given bipartite graph, and a 0*(2") time algorithm for computing the perma- 
nent of any n x n matrix. 

Theorem 2 Let F{xi, . . . , Xn) be any given Il„iY.sTl2 polynomial. It is #P-hard 
to compute the coefficient of any given multilinear monomial in the sum-product 
ofF. 

Proof It is well known (see Valiant [29]) that the problem of counting the 
number of perfect matchings in a bipartite graph is #P-hard. We shall reduce 
this counting problem to the problem of computing coefficient of a multilinear 
monomial in a polynomial. Let G = (Vi U V2, i?) be any given bipartite graph. 
We construct a polynomial F as follows. 

Assume that Vi = {wi, ■ ■ ■ , ft} and V2 = {ui, ■ ■ ■ , Ut}. Each vertex Wj S Vi 
is represented by a variable Xi, so is Ui S V2 by a variable yi. For every vertex 
Vi £ Vi, define a polynomial 

(vi,uj)eE 



Define a polynomial for the graph G as 

F{G)=F^---Ft. 

Let n = 2t, m = t, and s be maximum degree of the vertices in Vi. It is easy to 
see that F{G) is a rt-variate T\mX's^2 polynomial. 

Now, suppose that G has a perfect matching {xi,yi^), . . . ,{xt,yit). Then, 
we can choose ttj = Xjiji- from Fj, 1 < j <t. Thus, 

TT = TTi • 7r2 • ■ • TTt = a;ia;2 • • • Xtyiy2 ■■■yt 

is a multilinear monomial in FiG). Hence, the number of perfect matchings in 
G is at most c(7r), i.e., the coefficient of tt in F{G). On the other hand, suppose 
that F{G) has a multilinear monomial 



TT = TTj TTj = xia;2 • • ■ a;ij/iy2 • ■ • J/t 

in its sum-product expansion with tt' being a term from Fj, 1 < j < t. By 
the definition of _Fj, tt' = Xjyi-, meaning that vertices Vj and ui- are directly 
connected by the edge {j,ij)- Since tt' is multilinear, yi^, . . . ,yi^ are distinct. 
Hence, (xi, j/ij ),..., (a;t, j/i^) constitute a perfect matching in G. Hence, the 
coefficient c{'k) of tt in F{G) is at most the number of perfect matchings in G. 
Putting the above analysis together, we have that G has a perfect matching iff 
F(G) has a copy of the multilinear monomial tt = a;iX2 • • • Xtyiy2 ■ ■ -yt in its 
sum-product expansion. Moreover, G has c(7r) > many perfect matchings iff 
the multilinear monomial tt has a coefficient c(7r) in the expansion. Therefore, 
by Valiant's #P-hardness of counting the number of perfect matchings in a 
bipartite graph [5S], computing the coefficient of tt in F{G) is #P-hard. 



Theorem 3 There is a 0*(s(7i)3") time algorithm to compute the coefficients 
of all multilinear monomials in a polynomial F{xi, . . . ,Xn) represented by an 
arithmetic circuit G of size s{n). 



Proof We consider evaluating F from G via a bottom-up process. Notice that 
at most 2" many multilinear monomials can be formed with n variables. For 
each addition gate g in G with fan- ins fi, . . . , fs, we may assume that each fi is 
a sum of multilinear terms, i.e., products of distinct variables. This assumption 
is valid, because we can discard all the terms in fi that are not multilinear since 
we are only interested in multilinear monomials in the sum-product expansion 
of F. We simply add fi + ■ ■ ■ + fs via adding the coefficients of the same terms 
together. Since there are at most 2" many multilinear monomials (or terms), 
this takes 0(n2") times. 

Now we consider a multiplication gate g' in G with fan- ins hi and /i2. As 
for the addition gates, we may assume that hi is a sum of multilinear terms, 
i ~ 1,2. For each term tt with degree £ in /ii, we only need to multiply it with 



terms in /12 whose degrees are at most n — £. If the multiplication yields a non- 
multilinear term then that term is discarded, because we are only interested in 
multilinear terms in the expansion of F. This means that a term n of degree i 
in hi can be multiplied with at most 2"~^ possible terms in /12. Let rrii denote 
the number of terms in hi with degree i, 1 < i < n. Then, evaluating hi ■ /12 for 
the multiplication gate g' takes time at most 

0(n (mi 2"-i + 1712 2"-2 + . , . + rn^-i 2^)). (1) 

Since there are at most (") terms with degree i with respect to n variables, 
expression ((T|) is at most 

0{n [{^) 2"-i + (^) 2"-2 + , . . + (;j_ J 2"-"]) 

n 

= 0{n ^(r)2"-'') - 0{n 3"). 

Since C has s(n) gates, the total time for the entire evaluation of F for finding 
all its multilinear monomials with coefficients is 0(ns(n)3") = 0*(s(n)3"). □ 

The time bound in Theorem |3] can be improved when IISII polynomials are 
considered. 

Theorem 4 Let F{xi, . . . , Xn) be any given IlmT,sIlt polynomial. One can find 
coefficients of all the multilinear monomials in the sum-product expansion of F 
mO*(2") time. 



Proof Let F{xi, . . . ,Xn) = 111=1 -^» such that Fi = X) j=i "^jj ^^'^ "^ij i^ ^ 
term of degree at most t. We first consider Fm^i ■ F„i. Like what is done for 
the multiplication gate in the proof of Theorem [31 we multiply each term in 
Fjn-i with every term in Fj^. We discard all the resulting terms that are non- 
multilinear, because we are only interested in multilinear terms in F. Let G,„_i 
be the sum of all the remaining multilinear terms from Fm-i ■ F^. Then, G,„_i 
can have at most s^ < 2" many terms. Also, the time needed to obtain G,„_i is 
0{ts'^) = 0(is2"). Next, following the same approach, we do Fm-2 ■ Gm-i and 
let Gm~2 be the sum of all the remaining multilinear terms. The time needed 
to obtain Gm-2 is 0{ts2"). Continue this process to Fi ■ G2, we will have Gi as 
the sum of all the remaining multilinear terms that constitute all the multilinear 
monomials along with their respective coefficients in the sum-product expansion 
of F. The time for this last step also 0(ts2"). The total time for the entire 
process is 0(mis2") = 0*(2"). 



Corollary 5 There is a O* (1.415") tim,e algorithm to compute the exact num- 
ber of perfect matchings in a bipartite graph G = {Vi U V2,F) with n = 2\Vi\ = 
2IV2I vertices. 



Proof Let m = n/2, Vi = {«!,...,«„} and V2 = {ui,...,Um}. For each 
vertex Ui G V2, we define a variable Xi. For each vertex Vi E Vi, construct a 
polynomial 

where (vi,Ui.) G i? for j = I,- ■ ■ ,ii and Vi has exactly it adjacent vertices in 
G. Define 

H{G) = Hi---H,,/2. 

Then, H{G) is a (^)-variate n„/25^sni polynomial, where s = max{ii} < n/2. 
Following a similar analysis as in the proof of Theorem [5J G has a perfect 
matching iff if (G) has the multilinear monomial a;iX2 ■ • • a;„/2 in its sum-product 
expansion. Moreover, when there is a perfect matching, the number of perfect 
matchings in G is the same as the coefficient of a;ia;2 • • ■a;„/2. Therefore, by 
Theorem IH one can find the exact number of perfect matchings in G in time 
0*(2"/2) = C'*(1.415"). D 

The upper bound in Corollary[5]matches the best known deterministic upper 
bound of Ryser [26| for counting perfect matchings in a bipartite graph. The 
best known deterministic algorithm to compute the permanent of an n x n matrix 
is Ryser Algorithm [26] with 0*(2") time complexity that was devised almost 
50 years ago. A corollary of Theorem |4] implies an algorithm for computing the 
permanent of any matrix with the same time bound as Ryser algorithm does. 
Notice that when defining XISII polynomials in Section 2, we let the coefficients 
of all the terms in each clause to be 1 for ease of description. In fact. Theorems 
[3] and m still hold when arbitrary coefficients are allowed for terms in clauses of 
the input polynomial. 

Corollary 6 permanent The permanent of any given nxn matrix is computable 
in time O* {2"). 

Proof Let A = {aij)nxn be an n x n matrix with nonnegative entries Oij, 
^ l£ i,j ^ n. Design a variable Xi for row i and define polynomials in the 
following: 

Ri = {aiiXi-\ \-ainXn), 

P{A) = Ri---Rn- 

Let perm(A) denote the permanent of A. It follows from the above definitions 
that the coefficient of the multilinear monomial n = xi ■ ■ ■ Xn is precisely c(7r) = 
perm(A). Since R{A) is a n„I]„ni polynomial, by Theorem |4l we have the 
0*(2") time bound for computing perm(A). D 

The reduction in the proof of Corollary [5] implies the following result that 
somehow strengthens Theorem [21 

Corollary 7 It is #P-hard to computing the coefficient of any given multilinear 
monomial in an n-variate Hm^s polynomial. 
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4 Fully Polynomial-Time Approximation Schemes 
for HE Polynomials 

In this section, wc show that in contrast to Theorem [2] and Corollary [71 fully 
polynomial-time randomized approximation schemes ("FPRAS") exist for solv- 
ing the problem of finding coefficients of multilinear monomials in a nS poly- 
nomial and some variants of this problem as well. An FPRAS ^ is a ran- 
domized algorithm, when given any n-variate polynomial F and a monomial n 
together with an accuracy parameter e S (0,1], outputs a value A{F,'!T,e) in 
time poly(n, 1/e) such that with high probability 

(1 - e)c(7r) < A{F, IT, e) < (1 + e)c(7r). 

Theorem 8 There is an FPRAS for finding the coefficient of any given multi- 
linear monomial in a IlmT's polynomial F{xi, . . . ,.t„). 

Proof Let F{xi, . . . , a;„) = Y\T=i ^i ^^^^ that Fi = Y^^j=i ^ij with Si < s. No- 
tice that any monomial in the sum-product expansion of F will have exactly one 
variable from each clause Fi. This allows us to focus on multilinear monomials 
with exactly m variables. Let tt = Xi-^ ■ ■ ■ Xi^ be such a multilinear monomial. 
We consider how to test whether w is in F, and if so, how to find its coefficient 
c(7r). 

For each Fi, we eliminate all the variables that are not included in tt and 
let F/ be the resulting clause and F' = F[- ■ ■ F^ . If one clause F/ is empty, 
then we know that tt must not be a in the expansion of F' , nor in F. Now 
suppose that all clauses F^, 1 < i < m, are not empty. We shall reduce F' to 
a bipartite graph G = {Vi U V2,E) as follows. Define Vi = {vi, . . . ,Vm} and 
V2 = {iti, . . . , Um}- Here, each vertex Vi corresponds to the clause F/, and each 
vertex Uj corresponds to the variable Xj . Define an edge {vi , Uj ) in F if Xj is in 
F,. 

Suppose that tt is a multilinear monomial in F (hence in F'). Then, each 
Xi- in TT is in a distinct clause F4 . , 1 < j < m. This implies that edges {vt.,Ui), 
1 < j < m, constitute a perfect matching in G. On the other hand, if edges 
(vt , Ui .) , 1 < j < m form a perfect matching in G, then we have that Xi is in the 
clause Ft . Hence, tt = Xi^ ■ ■ ■ Xi^ is a multilinear monomial in F' (hence in F). 
This equivalence relation further implies that the number of perfect matchings 
in G is the same as the coefficient of the multilinear monomial tt in F. Thus, 
the theorem follows from any fully polynomial-time randomized approximation 
scheme for computing the number of perfect matchings in a bipartite graph, 
and such an algorithm can be found in Jerrum em at el. |17| . □ 

In the following we shall consider how to compute the sum S{F) of the 
coefficients of all the multilinear monomials in a HS polynomial F. 

Theorem 9 There is an FPRAS, when given any n-variate Tlm'^s polynomial 
F(xi, . . . ,Xn), computes S(F). 
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Proof Let F{xi,. . . ,Xn) = YVILiFi such that Fi = J2T=i^ij '^ith Si < s. 
Since every monomial in the sum-product expansion of F consists of exactly one 
variable from each clause Fj, ii m > n then F must not have any multilinear 
in its expansion. Thus, we may assume that m < n, because otherwise F will 
have no multilinear monomials. Let H = (xi + ■ • • + a;„). Define 

F'{xu...,xn) = ^•ff"-™ = Fi--.i^,„.i7"-™. 

Then, F' is a n„I]„ polynomial. For any given multilinear monomial 



in F with Xi- belonging to the clause Fj, 1 < j < m, let Xi^_^-^, . . . ,Xi^_^ be 
the n — m variables that are not included in tt, then 

/ 

^ "^^l ' ' ' '^im ' ^im + l ' ' ' '^in — i-n X\'^2 ' ' ' -^Tl 

is a multilinear monomial in F' . Because F' have n clauses with n variables, 
the only multilinear monomial that may be possibly contained in F' is the mul- 
tilinear monomial n' = a;ia;2 ■ ■ ■ Xn- If F' indeed has the multilinear monomial 
tt' with Xi- in the clause Fj, 1 < j < m, then tt = cc^^ • • • Xi^ is a multilinear 
monomial in F. This relation between tt and tt' is also reflected by the relation 
between the coefficient c{'k) of tt in the expansion of F and the efficient c(7r') of 
tt' in the expansion of F' . Precisely, the coefficient c(7r) of tt in F implies that 
there are c(7r) copies of x^j • ■ • Xi^ for the choices of the first m variables in tt'. 
Each additional variable Xi- , m + 1 < j < n — m, is selected from one copy of 
the clause H. Since H = {xi + ■ ■ -Xn), there are (n — m)! ways to select these 
(n — m) variables from [n — m) copies of H in F' . Hence, tt contributes a value 
of c(7r)(n — m)! to the coefficient of w' in F'. Adding the contributions of all the 
multilinear monomials in F to tt' in F' together, we have that the coefficient of 
TT in F' is S{F) ■ {n-m)\. By Theorem HI there is an FPRAS to compute the 
coefficient of tt' in F' . Dividing the output of that algorithm by [n — m)! gives 
the needed approximation to S{F). D 

We now extend Theorem |9] to HEn x HE polynomials. 

Theorem 10 Let F{xi, . . . ,Xn) be Ilk'SaTlt x ^m'^s polynomial with a > 2 
being a constant. There is a 0{a^ poly{n, 1/e)) time FPRAS that finds an e- 
approximation for the coefficient of any given multilinear monomial n in the 
sum-product F if n is in F, or returns "no" otherwise. Here, < e < 1 is any 
given approximation factor. 



Proof Let F ^ Fi- F2 such that Fi is a Hk^Jlt polynomial and F2 is a TLmTis 
polynomial. We first expand Fi into its sum-product expansion. Since we 
are only interested in multilinear monomials, all those that are not multilinear 
will be discarded from the expansion. We still use Fi to denote the resulting 
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expansion. We will have at most a'^ multilinear monomials in Fi as expressed 
in the following 

Fi=^6.V'^, (2) 

where bi = c(V') is the coefficient of the multilinear monomial ipi in F. 

Given any multilinear monomial tt, we consider how to test whether it is in 
F and if so, how to find its coefficient c{tt). Assume that tt is a multilinear 
monomial in F. Since F = Fi ■ F2, t: must be divided into two parts tt = tti • 7r2 
such that TTi is chosen from Fi and 7r2 is chosen from i^2- By expression ^, tx\ 
must be fAi for some \ <ij < a^ . If this not true, then tt is not in _F, so return 
"no". Now, for each -0^, such that -0^, is a possible candidate for tti, we decide 
whether 1x2 is a multilinear monomial in F2 and if so, we let 'K2{ipi) denote the 
second part of tt with respect to the first part tti = ipi- and find its coefficient 
c{TT2{ipi-)) in F2. By Theorem IH there is an FPRAS A to accomplish this task, 
since F2 is a lim^s polynomial. Let Aiipi) denote the approximation to the 
coefficient c(7r2(7ri.)) returned by the algorithm A with respect to the candidate 
ipi . . Let "tpi-^ , ■ • ■ , 'ipie be the list of all the candidates for tti . Then, the algorithm 
A returns A{tt) as 

^(tt) ^ b,,AiiJ^,) + ■ ■ ■ + b,,A{^p^,). 

Since A is an FPRAS, we have 

A{tt) < b,,{l + e)c(V'.,i • 7r2(V',;i)) + --- + b,,{l + e)c(0,;, • TT2{tl^,J) 

= (1 + e)[b^^c{'^pi^ •7r2(^ii)) H \- h^ci'ipi^ •7r2(V'iJ)] 

= (l + e)c(7r). 

Similarly, we have 

A{7r) > 611(1 - e)c(V'ii •7r2(V^iJ) H ^ b^^{l - e)c{ipi^ ■ Tr2{ipit)) 

= (1 -e)[&jic(0ii •7r2(V'ii)) H \-bii,c{i;i^ •7r2(V'jJ)] 

= (l-e)c(7r). 

Thus, ^(tt) is an e-approximation to c(7r). The time for expanding Fi is 
0{ta'^) ~ 0{na^). The time of the algorithm ^, by Thcorem[8l is 0(poly(n, 1/e)). 
So, the total time of the entire process is 0(a'^poly(n, 1/e)). □ 



5 Inapproximability 

Although in the previous section we have proved that there exist fully polynomial- 
time randomized approximation schemes for the problem of computing coeffi- 
cients of multilinear monomials in HmSs polynomials, yet in this section we 
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shall show that this problem is not approxiniable at all in polynomial time for 
IlmSsIIt polynomials with i > 2, unless P=NP. Thus, a clear inapproximability 
boundary arises between t = 1 and t = 2 for Ilm'^s^t polynomials. 

We consider a relaxed setting of approximation in comparison with the e- 
approximation in the previous section. Given any n-variatc polynomial F and 
a monomial tt together with an approximation factor 7 > 1, we say that an 
algorithm A approximates the coefficient ciir) in F within an approximation 
factor 7, if it outputs a value A{F, tt) such that 

- c(7r) <A{F,Tr) <7c(7r). 

7 

We may also refer ^ as a 7-approximation to c(7r). 

Theorem 11 No matter what approximation factor 7 > 1 is used, there is 
no polynomial time approximation algorithm for the problem of computing the 
coefficient of any given multilinear monomial in the sum-product expansion of 
a nmE3n2 polynomial, unless P=NP. 



Proof Let F{xi, . . . ,a;„) = Y\"^i Fi be a nml]3n2 polynomial. With loss of 
generality, we may assume that every term Tij in each clause F^ is a product 
of two variables. (Otherwise, we can always pad new variables to any given 
n^ 5^3112 polynomial to meet the above clean format.) It follows from Chen 
and Fu [5| that the problem of testing multilinear monomials in this type of 
polynomials is NP-complete. 

Let TT be any given multilinear monomial. Obviously, tt is in F iff its coeffi- 
cient c(7r) in F is bigger than 0. Thus, testing whether tt is in i^ is equivalent 
to determine whether the coefficient of tt in F is bigger than 0. 

Since every monomial in the expansion of i^ is a product of exactly one 
term from each clause Fi, all monomials in F must have the same degree 2m. 
If 2m > n, then there is no multilinear monomials in F. So we only need to 
consider the case of 2m < n. Let H ~ (xi + X2 + ■ ■ ■ Xn) and define 

F' = Fi • F2 • ij("-2™) (3) 

Then, the only multilinear monomial that F' may possibly have is ip — xix ^ 
2 ■ ■ ■ x„. If TT is a multilinear monomial in F with the coefficient c{tt) > 0, then 
following a similar analysis as wc did in the proof of Theorem |9] we have that tt 
contributes c(7r)(7i — 2m)! to the coefficient c(?/') of -0 in F' . This further implies 
that F has a multilinear monomial iff F' has the only multilinear monomial ip 
with its coefficient c{ip) = S{F){n — 2m)!. In other words, F has a multilinear 
monomial iff c(V') > in F' . 

Assume that there is a polynomial time approximation algorithm A to com- 
pute, within an approximation factor of 7 > 1, the coefficient of any given 
multilinear monomial in a Ilm'^3^2 polynomial. Apply A to F' for the multi- 
linear monomial ip. Let A{ip) be the coefficient returned by A for ip. Then, we 
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have 

- c(V') < A{ij) < 7 c(V'). 

7 

This means that F have a multihnear monomial iff A^ip) > 0- Hence, we have a 
polynomial time algorithm for testing whether F has any multilinear monomial 
via running A on ip in F' . However, this is impossible unless P=NP, because 
it has been proved in Chen and Fu ^ that the multilinear monomial testing 
problem for F is NP-complctc. D 

By Theorem [31 there is a fully polynomial-time randomized approximation 
scheme for the problem of computing the sum of the coefficients of all the multi- 
linear monomials in a H^E^ polynomial. However, when HnjS^Hj polynomials 
are concerned, even if s = 3 and t = 2, this problem becomes inapproximable 
at all regardless of the approximation factor. 

Theorem 12 Assuming P ^ NP, given any n-variate nmS3H2 polynomial F 
and any approximation factor 7 > 1 , there is no polynomial time approximation 
algorithm for computing within a factor of ^ the sum S{F) of the coefficients of 
all the multilinear monomials in the sum-product expansion of F. 



Proof Consider the same n-variate HmS3H2 polynomial F{xi, X2, ■ ■ ■ , Xn) as 
in the proof of Theorem [TlJ Define F' as in expression ^ . With a similar 
analysis, we have that F has multilinear monomials iff the coefficient of the 
multilinear monomial tp = X1X2 ■ ■ ■ Xn has the coefficient S{F) [n — 2m)l. That 
is, F has multilinear monomials iff the coefficient 0(7/') of ijj is bigger than zero in 
F' . Hence, like the analysis for Theorem llli any polynomial time approximation 
algorithm for computing the coefficient c(^) in F' can be naturally adopted as 
a polynomial time algorithm for the multilinear monomial testing problem for 
HmS3H2 polynomials. Since the latter problem is NP-complete (see Chen and 
Fu [5]), the former algorithm does not exists unless P = NP. 



6 Weak Inapproximability 

In this section, we shall relax the 7-approximation further in a much weak set- 
ting. Here, we allow the computed value to be within a factor of the targeted 
value along with some additive adjustment. Weak approximation has been first 
considered in our previous work on approximating the exemplar breakpoint dis- 
tance |10j and the exemplar conserved interval distance [7] between two genomes. 
Assuming P ^ NP, it has been shown that the first problem does not admit 
any factor approximation along with a linear additive adjustment [TU], while 
the latter has no approximation within any factor along with a 0{n^'^) additive 
adjustment [7j. Wc shall strengthen the inapproximability results of Theorems 
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[TT] and [T^ to weak inapproximability for computing the coefficient of any given 
multilinear monomial in a IIEn polynomials. But first let us define the weak 
approximation. 

Definition 13 Let Z be the set of all nonnegative integers. Given four func- 
tions f (x) , h{x) , a{x) and /3{x) from Z to Z with a{x) > 1, we say that h{x) is 
a weak (a(x), /3(x))-approximation to f{x), if 

TTiaxjo,^^^^^^^-^! </i(a;) <a(x) /(a;)+/3(a;). (4) 

I otyx) J 

Theorem 14 Let a{x) > 1 and I3{x) be any two polynomial time computable 
functions from Z to Z. There is no polynomial time weak {a{x), 13 (x))- approximation 
algorithm for computing the coefficient of any given multilinear monomial in an 
n-variate IlmT'3Ti2 polynomial, unless P=NP. 



Proof Let F{xi, . . . ,Xn) = YViLi^i t>c a nml]3n2 polynomial. Like in the 
proof of Theorem lll[ we assume without loss of generality that every term in 
each clause Fi is a product of two variables. We further assume that 2m > n, 
because otherwise there are no multilinear monomials in F. 

Choose k such that fc! > 2a{n + k)f3{n + fc) + P{n + k). Notice that finding 
such a fc < 2ri is possible when n is large enough, because both a and /3 are 
polynomial time computable. Let H = {xi-\-X2+- ■ ■ x„) and G = (yi+j/2 + - ■ ■ yk) 
with yi being new variables. Define 

F' = F- iy"-2™ .Qk^p^...p^^. jjn~2m . jjk (-g^ 

It is easy to see from the above expression ([5]) that F has a multilinear monomial 
iff F' has one. Furthermore, the only multilinear monomial that F' can possibly 
have is -tp = xi ■ ■ ■ Xn ■ yi ■ ■ ■ yk- 

Now consider that F has a multilinear monomial tt with its coefficient c(7r) > 
0. Since the degree of n is 2m, let Xi-^, . . . ,Xi,^_.^,,^ be the variables that are 
not included in tt. Then, the concatenation of tt with each permutation of 
Xij , . . . , Xi„_2^ selected from H"~^"^ and each permutation oi yi, . . . ,yk chosen 
from G'^ will constitute a copy of the only multilinear monomial ip in F' . Thus, 
TT contributes c{TT){n — 2m)! fc! to the coefficient c(V') of -0 in i^'. When all the 
possible multilinear monomials in F are considered, the coefficient of c(0) in F' 
is S{F){n — 2m)lkl. If F' has a multilinear monomial, i.e., the only one ip, then 
F has at least one multilinear monomial. In this case, the above analysis also 
yields ci^j) = 5(F)(n - 2m)!fc! in F' . 

Assume that there is a polynomial time weak (a, /3)-approximation algorithm 
A to compute the coefficient of any given the multilinear monomial in a 11^ 23112 
polynomial. Apply A to F' for the multilinear monomial ip. Let A{')p) be the 
coefficient returned by A for "0. Then, by expression (|4]) we have 

A{tlj) < a{n + k) c{tp) + I3{n + k) 
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= a{n + k) S{F) {n-2m)\ k\ + (3{n + k), (6) 

c{ip) - f3{n + k) 



(7) 



a{n + k) 
S{F) (n-2m)! fc! - /3(n + fc) 



a{n + k) 

When F docs not have any multihnear monomials, then F' docs not either, 
implying S{F) = 0. In this case, by the relation ([6]), we have 

A{^) < p{n + k). (8) 

When F has multilinear monomials, then F' does as well. By the relation (O, 
we have 

> S{F){n-2m)\k\-p{n + k) 
^ ' ~ a{n + k) 

fc! - f3{n + fc) (2a(n + k)P{n + k) + /3(n + fc)) - /3(n + fc) 
~ a{n + fc) a{n + fc) 

= 2/3(n + fc). (9) 

Since there is a clear gap between (— cxd, /3(n + fc')] and (2/3(ri + fc), +00), inequal- 
ities ([H]) and (ini) provide us with a sure way to test whether F has a multilinear 
monomial or not: If A{il)) > 2j3{n + fc), then F has multilinear monomials. If 
A{'ip) < (3{n + fc) then F does not. Since A runs in polynomial time, I3{n + fc) is 
polynomial time computable and fc < 2n, this implies that one can test whether 
F has a multilinear monomial in polynomial time. Since it has been proved in 
Chen and Fu [5] that the problem of testing multilinear monomials a nmS3n2 
polynomial is NP-complete, such an algorithm A does not exist unless P=NP. 

Combining the analysis for proving Theorems [12] and [TH wc have the fol- 
lowing weak inapproximability for computing the sum of coefficients of all the 
multilinear monomials in a IIEII polynomial. 

Theorem 15 Let a{x) > 1 and P{x) be any two polynomial time computable 
functions from Z to Z . Assuming P ^ N P , there is no polynomial time weak 
{a{x), I3{x))- approximation algorithm for computing the sum S{F) of the co- 
efficients of all the multilinear monomials in the sum-product expansion of a 
I{rrJ1^^2 polynomial F . 

7 The Maximum Multilinear Problem and Its 
Approximation 

Given any IIEn polynomial F{xi^ . . . , a;„) = Fi ■ ■ ■ F„i, F may not have any 
multilinear monomial in its sum-product expansion. But even if this is the case, 
one can surely find a multilinear monomial by selecting terms from a proper 
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subset of the clauses in F, unless all the terms in F are not multilinear or F 
is simply empty. In this section, we consider the problem of finding the largest 
(or longest) multilinear monomials from subsets of the clauses in F. We shall 
investigate the complexity of approximating this problem. 

Definition 16 Let F{xi, . . . , a;„) ~ Fi ■ ■ ■ F„i be a Ilm^s^t polynomial. Define 
MAX-SIZE{F) as the maximum length of multilinear m,onomials tt ~ pii-^ ■ ■ ■ iTi^, 
with TT^. in Fi^, 1 < j < k and 1 < ii < ■ ■ ■ < ik- Let MAX-MLM{F) to he 
a multilinear monomial n such that \tt\ ~ MX-SLZE[F), and we call such a 
multilinear monomial as a MAX-multilinear monomial in F . 

The MAX-MLM problem for an n-variatc nSII polynomial F is to find 
MAX-MLM(F). Sometimes, wc also refer the MAX-MLM problem as the prob- 
lem of finding MAX-SIZE(F). Wc say that an algorithm A is an approximation 
scheme within a factor 7 > 1 for the MAX-MLM problem if, when given any 
nsn polynomial F, A outputs a multilinear monomial denoted as A[F) such 
that MAX-SIZE(i^) < 7|^(F)|. 

Theorem 17 Let X > 2 be a constant integer. Let F be any given n-variate 
Ilm^s^x polynomial with s > 2. There is a polynomial time approximation 
algorithm that approximates the MAX-MLM problem for F within a factor of 
A. 



Proof Let F{xi, . . . , x„) = Fi ■ ■ ■ Fm such that each clause Fi has at most s 
terms with degrees at most A. Let M = Mi ■ M2 ■ ■ ■ Mk be a MAX-multilinear 
monomial in F. Without loss of generality, assume |Mi| > IM2I > ■■■\Mk\. 
We shall devise a simple greedy strategy to find a multilinear monomial tt to 
approximate AL . 

We first find the longest term tti from a clause Fi^ . Mark the clause Fi^ off 
in F. Let tt = tti. From all the unmarked clauses in F, find the longest term 
7r2 from a clause Fi^ such that 1^2 has no common variables in tt. Mark i^^^ off 
and let tt = tti • 7r2. Repeat this process until no more terms can be found. At 
this point, we obtain a multilinear monomial tt = tti • 7r2 • • ■ tt^ . 

Notice that each term in F has at most A variables. Each tt^ may share 
certain common variables with some terms in M . If this is the case, then ttj will 
share common variables with at most A terms in M . This means that we can 
select at least I > [|-] terms for n. The greedy strategy implies that 
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|Mi| + -.- + |A4| ^ \M\ 
A A ■ 

Hence, 

MAX-SIZE(F) = \M\ < A |7r| . 

Therefore, The greedy strategy finds the monomial tt that approximates the 
MAX-multihnear monomial M within the factor A. D 



Theorem 18 Let F{xi, . . . ,Xn) be any given n-variate II„i^s^t polynomial. 
Unless P = NP, there can be no polynomial time algorithm that approximates 
MAX-MLM(F) within a factor of n^^'"''/'^, for any e > 0. 



Proof We shall reduce the maximum independent set problem to the MAX- 
MLM problem. Let G — {V, E) be any given indirected graph with V = 
{vi, . . . ,Vn}. For each edge (vi,Vj) S E, we design a variable Xij represent- 
ing this edge. For each vertex Vi S V, let d{vi) denote the number of edges 
connecting to it and define a term T{vi) as follows: 

^^"'^ 1 (n(..,.,)e^ ^..) • {li;=r'''"' y^.) , if dM <n-l. 
We now define a polynomial F{G) for the graph G as 

F{G) = (T(«i) + ...+TK))". 

From the above definitions we know that all terms T{vi), 1 < i < n, have the 
same length n — 1. The number of new variables added to define F{G) is at 
most n{n — 1). 

Suppose that G has an independent set of k vertices u^j^ , . . . , Wi^. . Then 
there is no edge to connect Vi- and Vi^ for 1 < j,£ < k and j ^ £. This 
means that terms T{yi) and T{vif) do not have any common variables, so tt = 
T{vi^) ■ ■ ■ T{vi^) is multilinear with length k{n~l). On the other hand, suppose 
that we can choose terms T{vti), • ■ • ) '^(^t/) such that tt' = T{vti) ■ ■ ■ Tivtj) is 
multilinear. Then, there are no edges connecting any two pairs of vertices vt 
and Vt^ for I < j,i < k and j ^ i. This further implies that vertices vt-^, . . . , Vf, 
form an independent set of size / in G. Notice that \it'\ = f{n — 1). 

It follows from the above analysis that G has a maximum independent set 
of size K, iff F(G) has a MAX-multilinear monomial of length lC{n — 1). Assume 
that for any e > 0, there is a polynomial time algorithm A to approximate the 
MAX-MLM problem within an approximation factor of ri^^~')/^. On the input 
polynomial F{G), we can use A to find a multilinear monomial A{F{G)) that 
satisfies 

/C(n-l) < [n + n{n-l)]'^'-'^^^ AiF{G))=n^~' AiF{G)). (10) 
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It follows from above ((TU]) that 

K < .-rW^. (11) 

n — 1 

By (jlip . we have a factor n^~'^ polynomial time approximation algorithm for the 
maximum independent set problem. By Zuckerman's inapproximability lower 
bound of n^^'^ [31] on the maximum independent set problem, this is impossible 
unless P=NP. D 

Hdstad |16j proved that there is no polynomial time algorithm to approx- 
imate the MAX-2-SAT problem within a factor of ||. By this result, we can 
derive the following inapproximability about the MAX-MLM problem for the 
rim ^2 112 ■ Notice that Chen and Fu proved |S] that testing multilinear mono- 
mials in a IIS2 n polynomial can be done in quadratic time. 

Theorem 19 Unless P=NP, there is no polynomial time algorithm, to approx- 
imate MAXM-MLM{F) within a factor 1.0476 for any given Yim'^2 112 P'^^V' 
nomial F. 



Proof We reduce the MAX-2-SAT problem to the MAX-MLM problem for 
rim 1^2 112 polynomials. Let _F = i^i A ■ • • A F™ be a 2SAT formula. Without 
loss of generality, we assume that every variable Xi in F appears at most three 
times, and if Xi appears three times, then Xi itself occurs twice and Xi once. (It 
is easy to see that a simple preprocessing procedure can transform any 2SAT 
formula to satisfy these properties.) The reduction is similar to, but with subtle 
differences from, the one that was used in [8| to reduce a 3SAT formula to a 
Dm E3 02 polynomial. 

If Xi (or Xi) appears only once in F then we replace it by j/ii2/i2- When Xi 
appears twice, then we do the following: If Xi (or Xi) occurs twice, then replace 
the first occurrence by yiiyi2 and the second by yisya- If both Xi and Xi occur, 
then replace both occurrences by yiiyi2- When Xi occurs three times with Xi 
appearing twice and Xi once, then replace the first Xi by yiiyi2 and the second 
by yisyn, and replace Xj by y^iyis- 

Let G = Gi • ■ • Gm be the polynomial resulted from the above replace- 
ment process. Here, Gi corresponds to Fi with boolean literals being replaced. 
Clearly, i^ is a 11^^2112 polynomial and every term in each clause has length 2. 
For each literal Xi in F ^ let t[xi) denote the replacement of new variables for Xj. 
For each term T in G, i^^(T) denotes the literal such that T is the replacement 
of new variables for it. From the definitions of the replacements, it is easy to 
see that the clauses Fi^, . . . ,Fi_. in F are satisfied by setting literals Xi- € Fi. 
true, 1 < i < s,]S t: = t{xi.^) ■ ■ ■ t{xi^) is multilinear with t{xi) being a term in 
Gi , 1 < J < s. This implies that the maximum number of the clauses in F can 
be satisfied by any true assignment is /C iff a MAX-multilinear monomial in G 
has length 2/C. 
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Now, assume that there is a polynomial time approximation algorithm A 
to find a MAX- multilinear monomial in G within a factor of 1.0476 Apply the 
algorithm ^ to G and let A{G) denote the MAX-multilincar monomial returned 
by A. We have 

22 
2/C < 1.0476 A{G) < — Aifi), 

- 21 2 

Thus, we have a polynomial time algorithm that approximates the MAX-2-SAT 
problem within a factor of ly- By Hastad's inapproximability lower bound on 
the MAX-2-SAT problem [16], this is not possible unless P=NP. 

Khot at el. |19| proved that assuming the Unique Games Conjecture, there is 
no polynomial time algorithm to approximate the MAX-2-SAT problem within 
a factor of -^^. Notice that ^-^ > 1.0604 > |f > 1.0476. This tighter lower 
bound and the analysis in the proof of Theorem [TOl implies the following tighter 
lower bound on the inapproximability of the MAX-MLM problem. 

Theorem 20 Assuming the Unique Games Conjecture, there is no polynomial 
time algorithm to approximate MAXM-MLM)[F) within a factor 1.0604 for any 
given Ylm ^2 112 polynomial F. 

Remark. When the MAX-MLM problem is considered for n„jl]2n2 polyno- 
mials, Theorem[T7]gives an upper bound of 2 on the approximability of this prob- 
lem, while a lower bound of 1.0476 is given by Theorem [T9l assuming P ^ NP, 
and a stronger 1.0604 lower bound is derived by Theorem [20] assuming the 
Unique Games Conjecture. There are two gaps between the upper bound and 
the respective lower bounds. It would be interesting to investigate how much 
these two gaps can be closed. 
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